中 国生 态 农业 学 报 2017 年 6 月 第 25 丛 第 6 期 i 
Chinese Journal of Eco-Agriculture, Jun. 2017, 25(6): 795-804 Sw 


DOT: 10.13930/j.cnki.cjea.170113 
孙 琦 ， 白 中 科 ， 曹 银 贵 . 基于 生态 风险 评价 的 采 煤 矿区 土地 损毁 与 复 垦 过 程 分 析 [J]， 中 国生 态 农业 学 报 ，2017，25(6): 
795-804 


Sun Q, Bai Z K, Cao Y G. Process analysis of land destruction and reclamation in coal mining areas using ecological risk as- 
sessment[J]. Chinese Journal of Eco-Agriculture, 2017, 25(6): 795—804 


基于 生态 风险 评价 的 采 煤 矿区 土地 损毁 与 
复 垦 过 程 分 析 


孙 琦 !， 自 中 科 2"， 曾 银 贵 ” 


(1. 中 国 地 质 大 学 土地 科学 技术 学 院 ”北京 ”100083; 2. 国土 资源 部 土地 整治 重点 实验 室 ” 北 京 ”100035) 


摘 Wo nl We 
价 模型 分 析 采 煤矿 区 的 生态 风险 动态 变化 。 通 过 采用 最 小 损毁 累积 模型 测算 风险 源 累 积 0 利用 
感 影像 结 A rt eo ep pr eee i 
份 中 矿区 由 于 采矿 和 复姓 活动 带 来 的 生态 风险 变化 。 结 果 表 明 : 随 着 煤炭 产量 的 增加 , 2010 年 土地 损毁 的 累积 
影响 范围 比 2001 年 增长 了 7 095.17 hm?; 但 是 已 复 星 排 土 场 的 损毁 累积 危害 下 降 明 显 ， 并 且 随 着 矿区 扰动 区 
域 的 东 移 , 采矿 对 研究 区 西部 的 影响 也 在 减弱 。 经 过 10 年 的 土地 复明 与 生态 重建 措施 ,已 复明 的 排 土 场 生态 
风险 值 呈 下 降 趋势 ,生态 系统 趋 于 稳定 。 在 采矿 扰动 区 中 ,2001 年 处 于 中 等 风险 及 以 下 的 区 域 占 比 为 0.02%， 
到 2010 年 为 16.77%。 特大 型 采 煤 矿区 扰动 范围 较 大 ,但 是 采 后 土地 复明 有 助 于 降低 局 部 生态 风险 的 影响 。 通 
过 研究 矿区 土地 损毁 情况 、 复 垦 过 程 和 复明 后 的 状态 ， 分 析 不 同 区 域 生 态 风 险 的 动态 变化 ， 可 以 为 矿区 的 生态 
环境 治理 和 区 域 发 展 规划 提供 参考 ,为 矿区 管理 和 相关 决策 提供 科学 依据 。 
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Abstract: With the development of coal mining and land reclamation, the landforms and environmental conditions in mining 
areas have been disturbed. It is necessary to take positive initiatives to prevent eco-risks during mining exploitations and con- 
struct a key ecological civilization and socio-economic development. This paper used Pingshuo Coal Mine in Shanxi Province 
as a case study to analyze and build an evaluation model of dynamic ecological risks in coal mining areas. Using 30 m x 30 m 
grid units, the study evaluated, mapped and classified at unit grid scale ecological risks in the mining area. Based on the study, 
the ecological risk in the mining area varied with land destruction and reclamation processes. A model of the cumulative effect 
of land destruction was used to assess risk sources. Also remote-sensing images and field measured data were used to calculate 
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the values of ecological vulnerability index in the region. Through index calculations, the ecological risk values for the mining 
area were determined for 2001 and 2010. Also through the combination of existing remote-sensing techniques with field 
measurements, the values were reasonably up-scaled for the entire area with a highly visualized evaluation results. Then the 
effects of potential risks were determined and optimized land use forms and mining technology identified using the dynamic 
changes in ecological risk. The results showed that with increasing coal production, the cumulative affecting area of land de- 
struction increased by 7 095.17 hm? from 2001 to 2010. However, the cumulative effect decreased obviously for reclaimed 
dumps. With the continuous eastward open-pit mining, the effect of mining on the west region of the study area weakened. 
After 10 years of land reclamation and ecological reconstruction, ecological risks in reclaimed dumps declined and the eco- 
system seemly stabilized. The medium and low risk area accounted for 0.02% of the disturbed area due to mining in 2001 and 
up to 16.77% in 2010. This implied that land reclamation reduced the effects of ecological risk in the region. The combined 
risk value for open-cast and underground mining was unstable, with a value of ecological risk as low as 0.45 or as high as 0.80. 
The risk of partial uneven subsidence and soil erosion clearly increased in frequency and intensity in the area. The study 
demonstrated that the analysis of ecological risk could reflect the processes of land destruction and reclamation via mining, 
changes in ecosystem structure and function, and the spatial heterogeneity in space-time fabric. Hence such analysis was a 
reference base for planning environment recovery and regional development. As a source of risk for the regional ecosystem, 
land reclamation in mining areas should focus on increasing the ecological functions of reclaimed areas. 


Keywords: Ecological risk assessment; Ecological vulnerability; Land destruction; Land reclamation; Coal mining area 


作为 环境 风险 管理 的 重要 组 成 部 分 及 科学 支 ” 脆弱 区 , 矿区 在 开采 的 过 程 中 不 可 避免 地 造成 大 面 
撑 ， 生 态 风 险 评价 是 当前 环境 管理 领域 中 的 重点 和 积 地 表 挖 损 和 塌陷 mI、 多 次 损毁 、 土 壤 生 态 破 坏 、 
难点 0 中。 由 以 往 的 研究 中 可 以 看 到 ,矿业 开采 与 能 。 植被 受 损 、 水 土 流失 、 土 地 和 水 污染 等 问题 I。 本 
源 生产 等 会 引发 一 系列 的 生态 环境 问题 , 威胁 着 区  ” 文 以 山西 省 平 朔 矿区 为 例 , 利用 2001 年 和 2010 年 
域 生态 环 境 的 稳定 与 可 持续 性 "5l。 学 者 们 对 矿区 的 。 的 Landsat TM 影像 数据 ,结合 矿区 地 表土 地 损毁 特 
生态 风险 评价 研究 多 为 重金 属 污染 风险 评价 和 基于 ” 征 及 开采 和 复 垦 两 类 活动 带 来 的 生态 环境 变化 , 分 
Fa ei ee 析 一 定时 期 内 矿区 生态 风险 变化 的 过 程 特征 。 


本 风 了 和 的 可能 性, 了 


采 对 土地 扰动 的 增强 ,对 矿区 生态 风险 的 研究 也 。 工 研究 区 本 
a a em 平 逆 矿区 位 于 山西 省 北部 翔 州 市 平 重 区 内 ， 属 


温带 半 干 旱季 风气 候 , 年 平均 气温 4.8~7.8 'C, 年 降 
区 生态 风险 研究 的 趋势 -131。 : | 
泪 量 ~ 草原 植 类 J 地 带 ; 
矿业 开采 使 矿区 土地 利用 类 型 处 于 快速 变化 的 。 下 站 428.2~449.0 mm， 属 村 二 原 人 被 基于。 也 市 性 
土壤 为 栗 钙 土 和 栗 宰 士 , 土壤 有 机 质 含量 低 、 结 构 
动态 过 程 , 矿区 的 土地 损毁 作为 风险 源 ,其 对 生态 。 、 疡 他 能 力 弱 ”是 又 后 河 流域 十 壤 伪 地 最 为 严重 
统 的 潜在 危害 是 一 定时 期 内 的 累积 风险 使 用 。 ee 
in es 和 全 ee 的 地 区 ,同时 本 地 区 冬 春季 节 风 大 风 多 ,地 表 干 燥 。 
静态 的 景观 态 风 险 评价 者 果 不 适宜 表征 矿区 RN 本 地 区 生态 系统 抗 逆 能 力 差 ， 属 黄土 高 原 典 型 的 生 
系统 的 组 织 与 更 新 ， 有 必要 从 可 持续 的 角度 定 量化 态 脆弱 区 。 平 六 矿区 属 于 大 型 露 井 联 采 矿区 P] 现 
生态 风险 评价 指标 ， 反映 矿区 生态 风险 的 动态 变 包含 3 7 个 大 型 露天 矿 、 3 7 个 井 工矿 处 于 不 同 的 开采 
js 而 学 者 对 于 矿区 的 生态 环境 评价 多 关注 在 运用 生 阶段 ， 矿区 范围 近 160 km2。 地 貌 为 黄土 低 山 丘 陵 ， 
命 周 期 理论 对 不 同 发 展期 的 矿业 城市 进行 评价 ““， ”海拔 高 度 1 300~1 400 me 自 2001 年 以 来 , 矿区 的 生 
且 遥 感 数据 多 用 来 判断 景观 动态 变化 特征 “ ”… 从 ”产品 现 出 持续 稳定 增长 的 趋势 ， 原煤 产量 由 2001 年 
时 空 变化 趋势 来 分 析 矿 区 的 生态 风险 变化 的 研究 尚 ”的 2 180 万 t 提 高 到 2010 年 的 10 388 万 t 增长 4.77 
鲜 见报 道 。 在 全 球 环境 形势 及 我 国 经 济 发 展 与 生态 。 倍 , 年 均 增 长 1 174 万 t。 截至 2010 年 矿区 已 复 妊 的 
保护 并 存 的 制度 下 ,研究 矿区 土地 损毁 情况 、 复 星 土地 面 和 只 为 1 084.87 hm2， 占 已 有 排 土 场面 积 的 
过 程 和 复 垦 后 的 状态 ,对 潜在 复姓 区域 进行 生态 风 ”33.45%, 土地 复姓 垦 类 型 主要 以 乔木 林地 和 灌木 林地 
险 评 价 ， 可 以 为 煤矿 区 的 生态 环境 治理 和 区 域 发 展 ” 为 主 。 


规划 提供 决策 参考 。 1.2 ”数据 来 源 
我 国 煤炭 资源 大 部 分 处 于 干旱 、 半 干旱 的 生态 本 文 以 2001 年 和 2010 年 的 TM 影像 为 主要 数 
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据 源 , 分 别 为 2001 年 8 月 20 日 和 2010 年 7 月 12 
日 的 Landsat 5 平台 数据 ， 以 ENVI 4.8 软件 为 操作 平 
台 ， 对 影像 进行 预 处 理 ， 主 要 包括 大 气 辐 射 校 正 、 几 
何 校正 、 影 像 裁剪 等 操作 ， 并 采用 人 工 神经 网 络 
(artificial neural neb) 分 类 与 目 视 解 译 相 结合 的 方法 ， 
获得 研究 区 2001 年 和 2010 年 的 土地 利用 数据 ; 结 
合 研究 区 特点 ， 将 研究 区 土地 划分 为 耕地 、 草 地 、 
林地 、 城 镇 用 地 、 农 村 居民 点 、 交 通 运输 用 地 、 露 
天 采 坑 、 排 土 场 剥 离 区 、 工 业 场 地 。 此 外 ,本文 参 
考 朔 州 市 国民 经 济 和 社会 发 展 统计 数据 、《 户 州 市 土 
地 利用 总 体 规划 (2006 一 2020 年 阔 、 朔 州 市 2010 年 
土地 利用 变更 数据 等 ， 对 遥感 影像 数据 作 补充 ， 确 
保 研 究 数 据 可 靠 ; 其 他 数据 资料 包括 矿山 实地 调研 
数据 成 果 、 矿 区 1 : 2 000 地 形 图 等 。 
2 研究 方法 
2.1 研究 思路 

矿区 土地 按 不 同 的 演变 阶段 可 分 为 3 种 类 型 : 
原 地 貌 、 扰 动 地 和 复 垦 地 。 从 土地 生态 系统 的 稳 
定性 来 看 ， 原 地 貌 处 于 原始 的 生态 系统 环境 中 ， 
具有 一 定 的 稳定 性 ; 扰动 土地 因 人 为 活动 而 处 于 
快速 动态 变化 中 ; 复 垦 地 是 在 露天 采 空 区 或 原 地 
貌 上 快速 堆 双 形成 的 人 工 重 塑 地 貌 上 进行 土壤 重 
构 与 植被 重建 ， 生 态 系统 处 于 再 造 过 程 中 ， 具 有 
一 定 的 波动 性 。 同时 ， 由 于 平 朔 矿区 是 典型 的 露 井 
联 采 矿区 ， 开采 类 型 为 镶 主 型 和 印加 型 ， 采矿 作 
业 相 互 作 用 加 大 了 地 质 灾害 、 土 壤 侵 蚀 、 水 土 流 
失 等 风险 的 发 生 。 


风险 源 和 风险 受 体 是 生态 风险 评价 的 重要 组 成 ， 


在 矿区 生态 风险 因果 链 研究 的 基础 上 I， 从 风险 压 
力 的 释放 和 承受 风险 的 能 力 两 方面 构建 生态 风险 评 
价 指标 体系 , 即 土地 损毁 压力 和 生态 脆弱 度 。 将 风 
险 源 确定 为 矿区 土地 损毁 的 各 个 环节 ,土地 损毁 造 
成 的 环境 压力 越 大 ， 生 态 系 统 承 受 的 风险 性 就 越 
高 。 基 于 最 小 累积 阻力 的 连通 性 模型 近年 来 在 生态 
安全 评价 和 景观 格局 优化 等 方面 广泛 应 用 
可 以 定量 表征 矿区 生态 风险 源 与 受 体 间 的 风险 累积 
关系 。 

风险 受 体 在 生态 风险 评价 中 通常 被 认为 是 研究 
区 所 在 的 各 类 生态 系统 ， 可 以 通过 生态 脆弱 性 评价 
来 判断 生态 系统 在 特定 时 空 尺度 相对 于 外 界 干扰 所 
具有 的 敏感 反应 和 自我 恢复 能 力 王 ?。 如 果 生 态 环境 
本 身 较 脆弱 ， 其 抵抗 风险 和 恢复 的 能 力也 相应 较 低 ， 
承受 风险 的 能 力也 弱 ,， 选取 的 指标 既 要 反映 生态 系 
统 结构 的 健康 程度 ， 同 时 又 要 反映 生态 系统 服务 功 


能 的 强 弱 ,可 从 植被 、 气 候 条 件 、 土 壤 侵 刨 度 以 及 
环境 质量 等 方面 进行 评价 ?1。 
2.2 矿区 生态 风险 评价 方法 
2.2.1 评价 指标 体系 

基于 上 述 分 析 ,， 将 生态 风险 < 暴露 一 响应 ”评价 
过 程 具 体 化 为 对 风险 源 土地 损毁 累积 评价 和 风险 受 
体 生 态 脆 弱 度 评价 两 个 方面 ， 分 别 反映 “矿业 开采 - 
生态 系统 ”不 同 风险 源 、 风 险 受 体 及 其 交互 作用 的 特 
征 ， 可 以 表征 具体 的 开采 、 复 垦 活 动 或 生态 过 程 ， 通 
过 对 研究 区 30 mx30 m 的 像 元 进行 评价 实现 矿区 的 
生态 风险 空间 上 的 响应 。 

土地 损毁 累积 评价 的 准则 层 是 土地 损毁 压力 ， 
各 指标 通过 划分 矿业 活动 扰动 区 获得 矢量 数据 ， 进 
而 通过 空间 分 析 方 法 提取 指标 值 ; 二 级 指标 包括 土 
地 损毁 过 程 的 7 个 具体 指标 ( 表 1)。 风 险 受 体 生 态 脆 
弱 度 评价 则 包括 植被 、 湿 度 、 土 地 利用 类 型 等 5 个 
旨 标 ( 表 1)， 通 过 遥感 手段 和 矿区 野外 调研 实测 数据 
获得 。 

基于 德尔 菲 法 和 0~4 对 比 评分 法 确定 指标 权重 ， 
即 2 个 因子 比较 时 ,相对 重要 程度 有 以 下 3 种 基本 
情况 : 1) 影 响 很 大 的 因子 得 4 分 , 另 一 很 小 的 因子 得 
0; 2) 较 大 影响 的 因子 得 3 分 , 另 一 影响 较 小 的 因子 
得 1 分 ; 3) 同 样 影响 的 因子 各 得 2 分 。 对 这 12 个 指 
标 进 行 标准 化 , 定量 表征 矿区 生态 风险 ( 表 1)。 

对 指标 值 采用 极 差 标准 化 的 方法 ， 公 式 为 : 

正 效应 指标 : =(xo-xmin)/(xXmax—xmin) (1) 
负 效 应 指标 : XX,'=(xwmax-xa)/(xYmax-Xmin) ” (2) 

式 中 : 总 为 指标 xs 的 标准 化 值 , x 为 某 指标 x 的 取 值 ， 
xmax 为 指标 x 的 最 大 值 , x 为 指标 x 的 最 小 值 。 


2.2.2 ”评价 方法 
2.2.2.1 土地 损毁 累积 评价 


通过 土地 损毁 的 囚 积 作用 (cumulative effect， 

CE) 表 征 不 同类 型 损毁 土地 的 综合 危害 的 作用 大 小 ， 

通过 各 类 土地 损毁 压力 值 及 风险 源 对 外 扩散 耗费 系 

数 (sprawl resistance index, SRT) 进 行 量化 ， 即 风险 源 

作用 到 空间 上 某 点 的 生态 阻力 越 小 ， 则 此 点 累积 的 
土地 损毁 危害 越 高 ,生态 风 险 越 高 。 

CE; =1 一 2 (3) 


2 =V(1-p)xSRL (4) 
式 中 : CE, 表示 空间 内 任意 像 元 对 某 类 风险 源 的 累 
积 损毁 值 ，2' 表示 矿区 损毁 耗费 程度 的 标准 化 值 ， 
5 表示 对 应 像 元 的 土地 损毁 压力 指标 中 的 指标 权重 ， 
SRLI 是 该 类 风险 源 的 扩散 到 此 点 时 的 耗费 系数 的 标 
准 化 值 , i 为 30 mx30 m 的 像 元 单位 。 
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表 1 矿区 生态 风险 评价 指标 体系 
Table 1 Ecological risk evaluation index system of mining area 
准则 指标 定义 方向 ”权重 
Rule Index Definition Effect Weight 
土地 损毁 压力 ”地 下 开采 Underground mining 地表 和 裂缝 、 土 壤 退 化 Collapse-fissure, soil degradation 十 0.116 
了 P re of 下 
人 厂 石 排 痉 Gangue abandoned ”土壤 侵蚀 、 非 均匀 沉降 Soil erosion, unequal settlement 和 0.123 
录 | i 
ee 原 地 貌 剥 离 ee 原生 生态 系统 丧失 、 扬 侍 污 染 Loss of native ecosystems, dust pollution + 0.074 
Original landform stripping 
露天 开采 Open-pit mining 采 区 边 坡 稳定 性 、 粉 尘 污染 Instability of stripping slope, coal-dust pollution 十 0.098 
运输 排 弃 古 石 自燃 、 非 均匀 沉降 、 边 坡 稳定 性 、 土 壤 侵 蚀 、 扬 尘 污染 
D Gangue spontaneous combustion, unequal settlement, instability of mining + 0.291 
“pnB slope, soil degradation, dust pollution 
复 垦 管 护 Reclamation 植被 退化 、 土 壤 退 化 Vegetation degradation, soil degradation 证 0.091 
业 人 小 涝 涩 沁 沙 沁 沙 
洗 选 Washing 粉尘 污染 、 重金 属 污染 、 有 机 物 污染 0.208 
Pollution of dust, heavy metal and organic matter 
生态 脆弱 度 湿润 指数 Wet index 波段 差 值 提取 Calculation of different bands 一 0.202 
Ecological 
val rl 裸 土 指数 Bare soil index 波段 差 值 提取 Calculation of different bands 十 0.184 
index 植被 覆盖 度 Vegetation coverage 利用 NDVI 计 算 Calculation of NDVI 一 0.152 
生态 系统 服务 价值 生态 系统 服务 价值 当量 0 
Ecological system service value Ecosystem service value equivalent 
了 名 侧 ( 坡 度 、 植 被 类 型 量 娠 程度 
土壤 侵蚀 Soil erosion 水 蚀 + 风 蚀 (坡度 、 植被 类 型 、 年 降水 量 、 复 时 程度 ) | 二 
Analysis of slope, vegetation type, annual precipitation and extent of reclamation 
+: 正 效应 ; -: 负 效 应 。“+”means a positive effect, “—” means a negative effect. 
ne SRI 可 借助 GIS 中 Cost distance ” 土 像 元 , 植被 覆盖 度 (Pv) 采 用 下 式 计算 : 
a 抵抗 外 来 风险 。 供 鉴 景 和 中 ®) 
元 心 TTI ” 
NDVI 加 NDYImin 


地 利用 类 型 的 划分 情况 , 将 Cost 值 设置 为 林地 赋值 
1， 草 地 赋值 0.6， 耕 地 赋值 0.4， 居民 点 赋值 0.3， 裸 
地 赋值 0.2， 采 矿 用 地 与 交通 设施 用 地 赋值 为 0 中 。 
2.2.2.2 ”生态 风险 脆弱 度 评价 

风险 受 体 生态 脆弱 度 各 指标 的 选择 依据 和 计算 
方法 如 下 : 

1) 湿 度 指数 (WET): 通过 监测 矿区 地 面 土壤 湿 
度 及 其 变化 规律 ， 可 为 恢复 矿区 土地 原貌 、 恢 复 植 
被 覆盖 等 工作 的 展开 提供 基础 信息 卢 ]。 采 用 Crist 
提出 的 公式 提取 所 需 指数 I， 根 据 植被 指数 信息 和 
地 表 温 度 信息 , 采用 温度 植被 干燥 指数 法 反 演 得 到 
研究 区 的 土壤 湿度 信息 。 

2) 裸 土 指数 (NDSD: 由 纯 裸 土 指数 SI 和 建筑 指 
数 IBI 二 者 合成 ,是 由 于 本 研究 区 中 不 只 有 单纯 的 
裸 土 , 还 有 相当 一 部 分 的 工业 用 地 硬化 表面 , 具体 
指标 选用 及 计算 方法 参考 徐 泗 秋 提出 的 指数 计算 方 
法 R71。 

3) 植 被 覆盖 度 : 0 
植被 覆盖 度 可 以 作为 农业 、 环 境 穆 感 监测 和 评估 分 
析 的 定量 因 nn 
覆盖 度 更 能 反映 植物 的 功能 属性 和 生态 系统 的 健康 
程度 户 ]。 本 研究 借鉴 前 人 研究 结果 90， 将 NDVI>0.5 
的 像 元 认为 纯 植被 像 元 ， 当 NDVI<0.2 时 ， 作 为 纯 裸 


式 中 : NDVIns 取 值 为 0.5, NDVIwi, 取 值 为 0.2。 

4) 生 态 服务 价值 : 土地 利用 类 型 是 最 常见 的 景 
观 表 达 方 式 之 一 ， 可 以 直观 表达 人 类 活动 过 程 ， 借 
鉴 生态 系统 服务 价值 的 表述 方式 来 进行 赋值 。 有 学 者 
根据 地 域 特 征 提出 了 山西 省 的 生态 服务 功能 价值 的 当 
量 S 4, 而 矿区 生态 系统 具有 一 定 的 特殊 性 1， 通 
过 参考 前 人 研究 得 到 不 同 土地 利用 类 型 的 生态 系统 
服务 价值 ( 表 2)。 


表 2 平 朔 矿区 不 同 土地 利用 类 型 的 生态 系统 服务 价值 当量 


Table 2 Ecosystem service value of different land use types in 
Pingshuo Coal Mine 


土地 利用 类 型 生态 系统 服务 价值 当量 
Land use type Ecosystem service value (¥:hm 2.a !) 

林地 Forest land 2 416.00 

草地 Grassland 1 208.21 

耕地 Cultivated land 736.00 
裸 地 Unused land 0 
建设 用 地 Construction land 0 

采矿 用 地 Mining land —643.38 


5) 土 壤 侵 蚀 度 : 土壤 侵蚀 程度 反映 土 体 的 稳定 
性 。 由 于 采矿 、 排 土 造成 土壤 侵蚀 加 剧 , 研究 区 内 
土壤 侵蚀 包括 风蚀 、 水 蚀 和 工程 侵蚀 ,考虑 因子 包 
括 坡度 、 植 被 覆盖 度 等 。 根 据 平 朔 矿区 历年 统计 数 
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据 、 资 料 和 图 件 得 出 研究 区 内 不 同 区 域 的 土壤 侵蚀 
模 数 。 
生态 脆弱 度 的 计算 方法 如 下 : 
EVL= > WxDy (9 


式 中 : EVIL, 表示 露天 矿 内 任 一 像 元 上 的 生态 脆弱 度 
指数 (ecological vulnerability index)， 人 为 各 指标 
权重 ，V; 为 Vj 的 标准 值 , Vj 为 相应 的 指标 值 ,j 表示 
生态 脆弱 度 评价 中 的 评价 指标 , i 为 30 mx30 m 的 像 
元 单位 。 

2.2.3 矿区 生态 风险 评价 

借鉴 理想 景观 向 量 模型 54 可 得 到 生态 风险 评 
价 的 计算 方法 如 下 : 

ERV, = JEVI x CE, (7) 
式 中 : ERV; 是 研究 区 空间 单元 上 的 生态 风险 值 
(ecological risk value)，EVL 是 生态 脆弱 度 值 ，CE， 
是 土地 损毁 累积 作用 值 , i 为 30 mx30 m 的 像 元 单位 。 

为 了 进一步 说 明生 态 风险 值 的 空间 分 布 差异 
和 土地 扰动 类 型 间 的 变化 ,在 研究 区 设置 两 条 剖面 
线 ,经 过 原 地 貌 、 已 复 垦 排 土 场 、 工 业 场 地 、 未 复 
星 排 土 场 、 井 工 开采 、 露 井 复 域 、 露 天 采 区 和 剥离 
区 8 类 区 域 (图 1)， 对 不 同 区 域 的 生态 风险 值 进行 
比较 。 

3 结果 与 分 析 
3.1 矿区 土地 损毁 生态 风险 累积 

随 着 开采 活动 的 推进 ， 研究 区 土地 扰动 类 型 面 
积 扩大 、 类 型 增加 (图 1)。2001 一 2010 年 露天 开采 造 
成 大 量 土地 被 剥离 、 控 损 、 塌 陷 、 压 占 和 占用 ， 原 
地 貌 向 剥离 区 、 采 场 和 排 土 场 等 类 型 转变 ,新 增 土 
地 损毁 面积 为 2 767.79 hm”， 露天 采 场 面积 增加 
428.09 hm”, 未 复姓 排 土 场面 积 增加 为 1 281.32 hm”， 
其 中 二 次 损毁 的 排 土 场面 积 为 1 449.04 hm; 至 
2010 年 研究 区 已 复姓 排 土 场 1 084.87 hm2， 与 损毁 
土地 的 面积 比 为 1 : 4, 低 于 2001 年 的 1 : 3。 研 究 
期 内 年 均 复 垦 面 积 为 46.93 hm2， 年 均 土地 损毁 面积 
276.78 hm”， 复 垦 速 度 低 于 土地 的 损毁 速度 。 从 土地 
扰动 类 型 的 年 度 变 化 来 看 ,2001 一 2010 年 剥离 区 的 
变化 最 为 剧烈 ， 年 均 变 化 率 为 16.92%， 其 次 是 未 复 
坚 排 土 场 , 为 14.61%。 

2001 一 2010 年 土地 损毁 累积 危害 值 (CE) 的 空间 
分 布 、 区 域 大 小 发 生 了 明显 变化 (图 2)。2001 年 和 
2010 年 的 损毁 累积 危害 高 值 区 域 (0.7 达 CE<1) 为 剥 
离 区 、 采 矿区 和 工业 场地 ， 中 值 区 域 (0.35 达 CE<0.7) 
是 未 复明 的 排 土 场 和 距离 较 近 的 原 地 貌 ， 低 值 区 域 


(0 三 CE<0.35) 是 远离 矿业 开发 扰动 的 原 地 貌 。 随 着 
煤炭 开采 量 增加 及 井 工 煤矿 和 研究 区 东北 方向 东 露 
天 矿 的 投入 使 用 ，2010 年 土地 损毁 的 累积 影响 范围 
(CE>0.35) 达 24 113.23 hm”， 上 比 2001 年 增长 了 
7 095.17 hm>。 其 中 , 已 复姓 排 土 场 的 损毁 累积 值 明 
显 下 降 ,， 并 随 矿 区 扰动 区 域 的 东 移 ， 对 研究 区 西部 
的 影响 减弱 。 
3.2 ”矿区 生态 脆弱 度 指标 分 析 

将 矿区 生态 脆弱 度 评价 指标 分 为 5 个 等 级 ( 表 3)， 
定量 揭示 两 个 研究 期 5 个 指标 的 面积 变化 。2001 年 
WET 指数 主要 分 布 在 0.4~0.8，2010 年 则 主要 分 布 在 
0.6~1.0， 整 体 呈 增 长 趋势 ， 且 在 0.8~1.0 范围 内 2010 年 
比 2001 年 增加 8 394.91 hm”。NDSI 指数 在 2001 年 主 
要 分 布 在 0.4~1.0 范围 中 , 而 2010 年 向 低 值 区 域 转移 ， 
主要 分 布 在 0.2~0.8, 2010 年 0.4~0.6 中 值 面积 明显 增 
加 。 两 期 植被 覆盖 度 在 不 同 值 域 的 分 布 较 稳 定 , 2010 年 
相 比 2001 年 植被 覆盖 情况 总 体 呈 下 降 态 势 ， 植 被 覆盖 
度 小 于 20% 的 区 域 面 积 增长 较 多 。 矿 区 的 生态 服务 价 
值 低 值 区 面积 (0~0.4) 在 2010 年 相 比 2001 年 增长 了 1 
倍 ; 2001 一 2010 年 间 , 由 于 煤矿 的 开采 ， 原 地 貌 耕地 、 
林地 和 草地 面积 大 量 减 少 , 转化 为 剥离 区 、 露 天 采 坑 、 
排 土 场 等 土地 利用 类 型 ,生态 系统 服务 能 力 下 降 ; 
0.4~0.6 这 一 范围 的 面积 也 呈 下 降 趋势 ， 由 于 复明 地 多 
适宜 优先 恢复 为 林 草 地 ， 故 区 域内 耕地 面积 减少 ; 高 
值 区 数据 变化 不 明显 , 研究 区 内 林 草 地 的 损毁 与 恢复 
水 平 基本 持平 。 研 究 区 原 地 貌 是 风蚀 和 水 蚀 并 重 区 ， 
但 在 矿区 内 由 于 采矿 工程 造成 的 工程 侵蚀 重 于 原 地 
貌 ， 随 着 排 土 场面 积 的 扩大 ,土壤 侵蚀 强烈 的 区 域 面 
积 也 在 增加 (0.6~1)。 

结合 研究 区 的 扰动 变化 图 (图 1) 进 行 交 互 分 析 ， 
有 助 于 了 解 到 矿区 生态 环境 的 退化 和 改善 情况 。 从 
3 中 可 以 看 出 剥离 区 和 未 复 垦 排 土 场 等 区 域 其 
生态 脆弱 度 属于 高 值 区 ,但 是 会 随 着 矿业 开采 的 
推进 和 复 垦 工程 的 实施 而 发 生 空间 的 变化 。 这 两 类 
区 域 的 地 表 植 被 组 成 基本 没有 ,保水 能 力 差 , 土壤 
侵蚀 严重 ， 容 易 引 发 地 质 不 稳定 、 水 土 流失 、 扬 尘 
污染 等 问题 。 工业 场地 区 的 生态 脆弱 强度 没有 明显 
变化 ,在 空间 上 表现 为 面积 的 扩大 。 已 复 垦 的 排 土 
场 的 评价 结果 较 复 杂 ， 总 体 为 低 生 态 脆 弱 度 ， 并 表 
现 出 良性 变化 趋势 ， 局 域 生态 系统 稳定 性 优 于 原 
地 貌 。 其 中 已 复 垦 内 排 、 西 排 及 南 排 的 生态 系统 恢 
复 情况 极 好 ， 少 数 中 度 脆 弱 区 受到 早期 排 土 技术 
不 科学 和 后 期 管 护 不 当 的 影响 : 南 排 西 南 区 域 由 
于 古 石 自燃 造成 植被 大 量 死 亡 ， 大 面积 的 表土 裸 
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1 用 于 描述 2001 年 和 2010 年 矿区 扰动 区 域 生态 风险 变化 的 平 朔 矿区 土地 扰动 区 及 剖面 线 示意 图 
Fig.1 Sketch maps of mining disturbed zones and hatching lines in Pingshuo Mining Area used for analysis of ecological risk of 


disturbance in 2001 and 2010 
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土地 损毁 累积 
Land destruction 
cumulative effect 
mm 1 

“| 0 


2 2001 年 和 2010 年 平 朔 矿区 土地 损毁 风险 累积 危害 
Fig.2 Land destruction cumulative effect in Pingshuo Mining Area in 2001 and 2010 


露 , 还 有 部 分 区 域 有 灌浆 灭火 用 的 混凝土 浆 料 澄 
出 ; 西 排 部 分 区 域 由 于 火灾 的 影响 造成 了 植被 的 
退化 。 生态 脆弱 度 的 变化 反映 了 伴随 着 植被 消失 和 
重建 、 土 壤 结 构 发 生 较 大 改变 的 过 程 中 矿区 生态 系 
统 结构 及 功能 的 变化 。 
3.3 ”矿区 生态 风险 空间 异 质 性 

对 于 生态 风险 值 (ERV) 采 用 ArcGIS 中 的 自然 裂 


点 法 进行 划分 ，2001 年 的 裂 点 为 [0.00，0.13，0.26， 
0.40, 0.53，0.66, 0.79，1.00],2010 年 的 裂 点 为 [0.00， 
0.29, 0.40, 0.50, 0.61, 0.71, 0.82, 1.00]。 据 此 将 ERV 
分 为 4 个 等 级 : 高 风险 区 ,0.71<ERV 科 1; 较 高 风险 
区 , 0.50<ERV 夺 0.71; 中 等 风险 区 , 0.26<ERV 科 0.50; 
低 风 险 区 , 0.00 达 ERV 达 0.26)。 

从 图 4 可 以 看 出 ,2001 一 2010 年 平 朔 矿区 的 生 
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表 3 2001 年 和 2010 年 平 关 矿区 生态 脆弱 度 指 标 不 同等 级 的 面积 


Table 3 Areas of different levels of ecological vulnerability indexes in Pingshuo Mining Area in 2001 and 2010 hm” 


201711.02346Vv1 


chinaXiv 


指标 年 份 指标 值 范 围 Rang of index value 
Index Year 0~0.2 0.2~0.4 0.4~0.6 0.6~0.8 0.8~1 
湿度 指数 2001 0.90 23.60 14 644.52 17 011.14 18.66 
Wet index 2010 1.98 3.24 64.01 23 216.02 8 413.57 
裸 土 指数 2001 7.90 427.26 9 811.74 20 731.26 720.66 
Bare soil index 2010 28.16 1 599.04 19 580.32 10 486.04 5.26 
植被 覆盖 度 2001 23 905.38 3 599.30 1 889.00 1 154.55 1 150.59 
Vegetation coverage 2010 26 089.93 2 662.16 1 203.06 665.34 1 078.33 
生态 服务 价值 2001 768.61 2 645.44 16 062.31 8 268.56 3 957.27 
Ecological system service value 2010 1 588.78 4 912.80 12 584.30 8 708.82 3 904.22 
土壤 侵蚀 度 2001 11 812.90 17 719.70 664.30 30.76 1 471.79 
Soil erosion 2010 10 677.16 13 977.66 2 934.11 2 029.11 2 082.92 


生态 脆弱 度 
Ecological 
vulnerability 


2001 2010 


3 2001 年 和 2010 年 平 朔 矿区 生态 脆弱 度 空间 分 布 变 化 
Fig.3 Spatial distribution maps of ecological vulnerability in Pingshuo Mining Area in 2001 and 2010 
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Ecological risk 
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4 2001 年 和 2010 年 平 朔 矿区 生态 风险 分 布 
Fig.4 Distribution of ecological risk in Pingshuo Mining Area in 2001 and 2010 
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态 风 险 分 布 发 生 了 较 大 改变 。2001 年 受 采 矿 扰动 区 
域 (采矿 区 和 复 垦 区 ) 的 像 元 数 为 34 377 个 ，2010 年 
为 105 964 个 , 分 别 占 两 期 研究 区 像 元 总 量 的 9.76% 
和 30.08%。 在 扰动 区 中 , 2001 年 处 于 中 等 风险 及 以 
下 的 区 域 仅 为 0.02%，2010 年 为 16.77%, 可 以 看 出 
虽然 矿区 开采 范围 在 扩大 , 但 是 随 着 采 后 复 垦 工作 
的 开展 , 局 部 生态 风险 降低 。 在 两 个 评价 期 中 均 为 
已 复 垦 排 土 场 的 安 太 堡 南 排 土 场 和 西 排 土 场 ， 生 态 
风险 均值 发 生 显著 变化 ， 南 排 土 场 从 0.64 降低 0.50， 
西 排 土 场 从 0.62 降低 到 0.46。 排 土 场 的 均值 差异 反 
映 平 半 矿区 多 年 复姓 期 的 变化 ,由 于 早期 排 土工 艺 
不 规范 南 排 土 场 出 现 诸如 滑坡 、 古 石 自燃 等 生态 问 
题 ， 评 价 结果 表现 出 较 高 生态 风险 ; 西 排 土 场 排 土 
和 复 盟 时间 晚 于 南 排 土 场 , 但 是 排 土 工艺 较 南 排 土 


1.0 
0.8 
0.6 
0.4 


0.2 


生态 风险 值 Ecological risk value 


加 ji : : 
0 50 100 150 200 250 


场 科学 合理 ， 对 生态 系统 的 维护 有 一 定 正 向 影响 ， 
经 过 10 年 的 管 护 表 现 为 低 生 态 风 险 。 

剖面 线 a 和 b 分 别 经 过 原 地 貌 、 已 复姓 排 土 场 、 工 
业 场 地 、 未 复 星 排 土 场 、 井 工 开采 、 露 井 复 域 、 露 天 采 
区 和 剥离 区 8 类 扰动 区 (图 1)， 其 经 过 的 像 元 风险 值 形 
成 剖面 曲线 ， 见 图 5。 剖 面 线 a 经 过 的 生态 风险 值 由 图 
5A 和 图 SC 表示, 图 中 显示 原 地 貌 的 生态 风险 基本 无 变 
化 ; 已 复姓 的 南 排 在 2010 年 风险 值 上 升 且 上 下 波动 ， 
这 是 由 于 南 排 发 生 了 大 面积 的 自燃 ， 局 部 生态 受 损 严 
重 ; 2010 年 已 复姓 的 安 太 堡 内 排 生 态 风险 值 明 显 下 降 ; 
2001 年 尚 在 进行 开采 和 剥离 的 区 域 已 经 发 展 成 为 工业 
场地 和 未 复 晨 的 排 土 场 ，2010 年 的 生态 风险 值 较 十 年 
前 有 小 幅 下 降 ; 未 复 晨 排 土 场 的 风险 值 远大 于 已 复 星 
的 排 土 场 。 


0 50 


像 元 顺序 Pixel order 


5 2001 年 (A, B) 和 2010 年 (C, D) 平 关 矿 区 生态 风险 剖面 a(A, C) 和 剖面 b(B, D) 分 析 
Fig.5 Ecological risk analysis along hatching line a (A, C) and hatching line b (B, D) in Pingshuo mining area in 2001 (A, B) and 
2010 (C, D) 
0O: 原 地 貌 ; R: 已 复姓 排 土 场 ; U: 未 复 垦 排 土 场 ; I: 工业 场地 ; OM: 露天 采 区 ; UM: 井 工 采 区 ; S: 剥离 区 ; OUM: 露 井 复 域 。 〇 D 安 太 堡 南 
排 ; @O 安 太 堡 内 排 ; @ 安 太 堡 西 排 ; 安家 岭 西 排 ; 安家 岭 东 排 。O: original land; R: reclaimed dump; U: un-reclaimed dump; I: industrial site; 
OM: open-cast mining area; UM: underground mining area; S: stripping area; OUM: open and underground combined mining area. () south dump 


of Antaibao; @) inner dump of Antaibao; @) west dump of Antaibao; © west dump of Anjialing; @ east dump of Anjialing. 


剖面 线 b 在 2001 年 所 经 过 的 区 域 相对 简单 ， 
2010 年 此 区 域 的 扰动 类 型 随 着 采矿 推进 而 丰富 起 来 
(图 SB 和 图 5D)。 位 于 研究 区 西南 部 的 原 地 貌 区 域 
在 2010 年 土地 利用 类 型 没有 发 生变 化 ,但 受到 井 工 


开采 的 影响 ,其 生态 风险 值 上 升 ， 且 不 稳定 ， 在 0.65 
的 数值 上 下 波动 ; 工业 场地 的 风险 值 基本 平稳 ; 安 
家 上 岭 西 排 经 过 10 年 的 土地 复 垦 措施， 生态 风险 值 显 
著 降 低 , 但 风险 值 变化 不 稳定 且 浮 动 较 大 , 低 值 小 
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(不甘 局 03 


于 0.45 而 高 值 达到 0.80， 该 区 域 受 井 工 开采 影响 属 
于 露 井 复 域 ， 增 加 了 如 排 土 场 的 不 均匀 沉降 、 边 坡 
稳定 性 、 土 壤 侵 刨 度 等 风险 发 生 的 概率 和 强度 ; 已 
复 垦 的 安家 岭 东 排 由 于 复 垦 年 限 较 短 且 距 离 露 天 采 
区 近 ， 生 态 风险 值 仍 较 高 ; 2010 年 处 于 两 处 剥离 区 
之 间 的 原 地 貌 生 态 风 险 值 受到 距离 剥离 区 的 远近 影 
响 而 出 现 高 低 不 同 的 变化 。 

4 结论 与 讨论 

1) 本 研究 以 位 于 生态 脆弱 区 的 平 朔 矿区 为 研究 
对 象 , 通过 构建 生态 风险 评价 模型 表征 煤矿 开采 过 
程 中 潜在 生态 风险 问题 以 及 土地 复 垦 与 生态 重建 工 
程 对 区 域 生 态 环境 的 保护 和 改善 。 该 模型 从 土地 损 
毁 累 积 危 害 和 生态 系统 脆弱 度 两 方面 评价 了 矿区 生 
态 系统 受 损 度 和 承受 风险 的 能 力 ， 综 合 反 映 矿区 经 
过 多 年 开采 和 复 垦 不 同 区 域 生 态 风险 的 变化 ,其 评 
价 结果 符合 矿区 的 实际 。 通 过 综合 遥感 技术 和 实地 
测量 数据 ,不 仅 可 以 定量 评价 矿区 生态 风险 值 ， 并 
实现 了 评价 结果 和 变化 的 可 视 化 。 

2) 矿 区 生态 风险 的 大 小 和 变化 与 矿区 开采 速 
度 、 土 地 损毁 面积 大 小 、 开 采 过 程 科学 性 和 土地 复 
妊 措施 密切 相关 。 由 评价 结果 可 知 ， 经 过 10 年 的 土 
地 复 垦 与 生态 重建 措施 , 已 复明 的 排 土 场 生 态 风险 
值 呈 下 降 趋 势 ， 生态 系统 趋 于 稳定 。 特 大 型 露 井 联 
采矿 区 扰动 范围 较 大 , 但 是 随 着 采 后 复 垦 工作 的 开 
展 ， 有 助 于 降低 局 部 生态 风险 的 影响 。 

3) 通 过 研究 矿区 的 土地 损毁 与 复 垦 过 程 , 评价 生 
态 风险 的 动态 变化 ,可 以 识别 和 确定 不 同 阶 段 的 潜在 
风险 影响 以 及 部 分 有 待 优化 的 用 地 类 型 和 采矿 环节 。 
在 开采 之 前 预测 可 能 发 生 的 生态 风险 , 在 采矿 过 程 中 
介入 人 工 干 扰 手段 尽量 避免 或 减少 开采 对 生态 系统 的 
破坏 ， 在 采矿 活动 结束 后 及 时 进行 生态 系统 的 恢复 与 
重建 , 可 以 有 效 避 免 生 态 系统 的 灾难 性 破坏 。 作 为 区 
域 生 态 系统 的 风险 源 , 矿区 土地 复 垦 应 从 单一 的 土地 
利用 类 型 的 数量 和 质量 需求 向 生态 功能 需求 转变 ， 如 
何 建 立 可 持续 发 展 的 生态 矿区 、 构 建 系统 性 和 统筹 规 
划 的 复 垦 矿区 有 待 进一步 的 研究 并 具有 现实 意义 。 
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